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The Nikolaevskiy model for pattern formation with continuous symmetry exhibits spatiotemporal 
chaos with strong scale separation. Extensive numerical investigations of the chaotic attractor reveal 
unexpected scaling behavior of the long-wave modes. Surprisingly, the computed amplitude and 
correlation time scalings are found to differ from the values obtained by asymptotically consistent 
multiple-scale analysis. However, when higher-order corrections are added to the leading-order 
theory of Matthews and Cox, the anomalous scaling is recovered. 
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I. INTRODUCTION 

The study of dynamics and pattern formation in non- 
linear spatially extended dynamical systems continues to 
attract considerable attention ^'-'s'l, with potential solu- 
tion behaviors ranging from attracting steady states to 
complex spatiotemporal dynamics including spatiotem- 
poral chaos (STC), which is characterized by temporally 
chaotic dynamics with decaying spatial correlations and 
a finite density of positive Lyapunov exponents [H, Q . 

In the search for general principles, much effort has 
gone into the derivation and study of simplified model 
equations to highlight and clarify pattern-forming prop- 
erties of more complicated "full" systems. Such paradig- 
matic models, such as the Ginzburg-Landau equation 
for the slow evolution of modulations of an underlying 
pattern, have been remarkably fruitful in predicting the 
leading-order scaling behavior and stability of patterns, 
by comparison with numerical or experimental results. 

The identification of relevant scaling laws forms an 
integral part of the asymptotic procedures used to de- 
rive such reduced model systems. The successful appli- 
cation of such methods rests on the expectation that the 
leading-order truncation obtained subject to a consistent 
scaling assumption indeed captures the dominant scaling 
of the full problem, in the sense that the perturbative 
effects of neglected terms are "small" and occur at the 
scaling implied by the dominant balance. In this paper, 
we present an example which appears to challenge this 
expectation. 

It is known from various examples that the stability 
of finite- wavelength pattern or wave phenomena may be 
complicated through interaction with a mean flow; for in- 
stance, in convection with free-slip boundary conditions 
[1| , in some reaction-diffusion systems @ and elsewhere, 
such coupling can lead to unusual instability mechanisms 
of the underlying pattern or flow, whose analysis may dis- 
play mixing of scales in the perturbation parameter and 
require higher-order terms in the perturbative expansions 
for a complete understanding Q. 



In the present work we study a model for one- 
dimensional pattern formation in the presence of contin- 
uous symmetry, for which the interaction of modes with 
two well-separated length scales leads to spatiotemporal 
chaos. Focussing in this case on characterizing the long- 
term chaotic regime rather than on its onset via the ini- 
tial stability problem, for the statistics of the long-wave 
mode we find a scaling discrepancy between predictions 
of leading-order multiple-scale analysis and computations 
of the full system. Unexpectedly, though, the agreement 
in scaling is restored by including next-order corrections 
in the amplitude equations. 

We investigate the Nikolaevskiy partial differential 
equation (PDE), given in canonical (derivative) form as 
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(where Ux = dxU = du/dx), which was originally pro- 
posed as a model for longitudinal seismic waves in vis- 
coelastic media Q, and has more recently been ob- 
tained in the context of phase dynamics in reaction- 
diffusion systems [13, [Hi and finite-wavelength trans- 
verse instabilities of traveling fronts jTl]. In fact, the 
Nikolaevskiy equation has a continuous symmetry — in 
the form ([T]) it is Galilean-invariant, and preserves the 
spatial mean — and appears to be a paradigmatic model 
for short-wave pattern formation with such symmetry 

Considering individual Fourier modes of the form 
exp(CTt -|- ikx), the linear dispersion relation about the 
trivial solution u = is cr(fc) = fc^[r — (1 — fc^)^]. Since 
for r < 0, all initial conditions of ([Ij decay to a spatially 
homogeneous solution, in the following we consider only 
r > and write r = e^; as usual the 0{e) width of the 
unstable band and 0{£'^) maximum growth rate suggest 
natural slow length and time scales X — ex, T — e^t. 
For < r = ^ 1, the Nikolaevskiy equation ([T]) has 
a one-parameter (up to translation) family of stationary 
periodic solutions, or rolls, Uq{x) for \q\ < 1/2, which 
may be found by weakly nonlinear analysis to be IT^ 

uq{x) = 6e^/T^^e'^^+'"^'^'' + c.c. + 0{e^). (2) 

Previous studies of dH) have revealed some surprises. 
Due to the presence of a symmetry-induced long-wave 
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(Goldstone) mode Uq interacting with the weakly unsta- 
ble pattern mode C/i , all roll equilibria Uq{x) are unstable 
for all r > pa - [l5j : instead, there is a direct supercrit- 
ical transition from spatial uniformity to spatiotemporal 
chaos [16]. Similar behavior occurs, for instance, in elec- 
troconvection in liquid crystals [l7|. Tsjl. and has been 
termed "soft- mode turbulence" (tI. [l5|. 

By deriving leading-order modulation equations using 
multiple-scale analysis, Matthews and Cox [l^ further- 
more showed that in the chaotic state, Ui (and hence u) 
has a rather unusual r'^/^ = £^/^ scaling. In the present 
work, through extensive numerical studies we show that 
Nikolaevskiy STC exhibits yet more unexpected scaling: 
amplitudes and correlation times of the large-scale mode 
Uq have an anomalous dependence on r, seemingly incon- 
sistent with leading-order multiple-scale analysis, though 
recovered with the addition of next-order terms to the 
amplitude equations. In particular, in the chaotic regime 
the mean-square amplitude of Uq appears to scale as r^/^ , 
rather than the predicted r^. 
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FIG. 1: (Color online) Time-averaged Fourier power spectrum 
S{k)/£ for the Nikolaevskiy equation ^ ior L = e£ = 27r-25.6 
and various £^ values. From top to bottom, the spectra shown 



are for e 



1.0, 0.1, 0.04, 0.01, 10" 



and 10" 



Note the 



separation of scales for s < 0.04 (blue solid curves) 



II. NUMERICAL METHODS, SOLUTION 
PROPERTIES AND MODAL DECOMPOSITION 

We numerically computed (mean-zero) solutions of the 
Nikolaevskiy equation ([Ij on an ^-periodic domain in the 
chaotic regime, over a wider e-range than in earlier stud- 
ies: 10~^ < £^ < 0.2. The system size £ was chosen 
so that the number of Fourier modes in the unstable 
band of the dispersion relation cr(fc), about [ef/27rj, suf- 
ficed to capture the large-system limit; to minimize sys- 
tem size effects across our computations, we fixed the 
domain length in the slow space variable X = ex, at 
L = ee = 2n ■ 128m/10 with m=l and 2 (about 12 and 
25 unstable modes, respectively), and checked our results 
with TO=4. We used a pseudo-spectral method in space 
with N = 2"^ Fourier modes for J=ll to 16, chosen so 
that the maximum wave number retained was well within 
the strongly decaying regime of cr(fc): — N ■2tt/£> 6. 
Our time integration was performed using a fourth-order 
Runge-Kutta exponential time differencing (ETDRK4) 
scheme [H,!!^], with a typical time step of At — 0.02/e^, 
fixed in the slow time variable T — e^t; after a transient, 
we computed statistics over a time interval of length 

tra = T,„/£2 = 4 • IOVe'- 

The (time-averaged) power spectrum S{k) = (|ufcp) of 
Fig. [1] captures the dominant features of the numerically 
observed chaotic dynamics: As pointed out by Tanaka 
[2ll |. for "large" e the STC appears qualitatively similar 
to that of the related well-known Kuramoto-Sivashinsky 
equation Fourier modes with < fc < 1 are (at most) 
weakly damped, and the power spectrum reveals no sig- 
nificant scale separation for > 0.1 (dashed curves). 

However, for sufficiently "small" e, we observe the dis- 
tinct gap between the fc « and « 1 energetic modes 
which characterizes "Nikolaevskiy chaos" ^2lJ . Note that 
for £ <C 1, the subdominant peaks at |fc| « 2,3,... repre- 




FIG. 2: Gray-scale visualization of a computed solution of the 
Nikolaevskiy equation ((T)) for r — = 0.04, with domain size 
I — 1600 and over a time interval of length t = 400; since the 
roll wavelength \ 2-k <^ £, the small-scale rolls are barely 
visible in this representation. 



sent higher harmonics slaved to the Fourier modes |fc| « 1 
via the quadratic nonlinearity; we remark that exponen- 
tial time differencing integration schemes [is| are partic- 
ularly well-suited to capturing this slaving behavior and 
improving the accuracy at small scales. 

The strong scale separation of the spatiotemporally 
chaotic dynamics in this small-e regime is readily ap- 
parent in a representation of a typical numerically com- 
puted time series on the multiple-scale chaotic attractor 
of the Nikolaevskiy PDE ([T|), as in Fig. [51 or of a snap- 
shot of a typical solution as seen in Fig. Efa): Solutions 
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FIG. 3: (Color online) (a) Snapshot at a fixed time ti of a 
solution u{x,ti) (blue) of the Nikolaevskiy equation ^ with 
r — — 0.04 and domain length £ — 800. The figure also 
shows the reconstructed mean mode Uo (black) and enve- 
lope Uo ± 2|J7i| (black dashed line), (b) Instantaneous spec- 
trum |Tifc(ti)P of the solution in (a), showing the wavenumber 
ranges for the Fourier filters. 



on the attractor are characterized by long-wave modula- 
tions, aperiodic in space and time, of small-scale rolls of 
wavelength X ^ 2tt(1 + 0{e)) (corresponding to the lin- 
early unstable modes of cr(k)). This motivates a natural 
decomposition of the solution u{x,t) of the Nikolaevskiy 
PDE dll): with X = ex,T = e^t, we have 

uix,t) = C/i(X,r)e" +C.C. + ;7o(X,T) + h.o.t., (3) 

where Ui and Uq are slowly varying envelopes of the 
dynamically significant Fourier modes near = 1 and 
k = 0, respectively. 

We have implemented Fourier filters to extract Uq and 
Ui at each time step from the Fourier space numerical 
solution Uk{t) of the full PDE (P; see Fig. Efb): For 
some fc* £ (0, 1), we identiiy the modes with < A:* as 
being "large-scale", while those with \k\ > fc, belong to 
the patterned mode and its higher harmonics. In view 
of the strong suppression of the linearly damped modes 
near k — 0.5 for < 0.04 (see Fig. [TJ, we may safely 
choose A;* — 0.5; we have verified that our results are 
insensitive to the choice of fc*. Thus we obtain Uq from 
{uk\ \k\ < 0.5}, and C/ie" and hence Ui from {uk\ 0.5 < 
k < 1.5}. In Fig.|31[a) we have also shown the mean mode 
Uq and the reconstructed envelope Uq ± 2\Ui\. 



of the finite-wavelength spatial pattern mode with the 
lon g-wa ve mode Uq induced by the continuous symmetry 
1311231 [23 | . In this context it is of particular interest to 
determine the exponents a and /3 describing the scaling 
of Ui and Uq as e ^ 0. Thus one postulates the Ansatz 



u{x, t) - e°'A{X, T)e" + c.c. + e^/(X, T) -f 



(4) 



where A and / are assumed 0(1), and / has mean zero. 

The 0{£) amplitude of the known roll solutions Uq{x) 
^ suggests the apparently natural choice a=l (see for 
instance familiar from other pattern- forming sys- 

tems [l|, which we will denote the Ginzburg-Landau (GL) 
scaling. It turns out, though, as discussed in [l^, that 
amplitude equations for ([1]) derived using a=l in ^ are 
asymptotically inconsistent for any choice of /3 if / ^ 0. 

Instead, as shown by Matthews and Cox ^SJ, the only 
asymptotically consistent amplitude equations are ob- 
tained by using the Ansatz @ with a=3/2, /3=2; we shall 
call this the MC scaling. Standard multiple-scale tech- 
niques then yield the leading-order amplitude equations 
proposed to describe STC in the Nikolaevskiy equation 
dl]): the Matthews-Cox (MC) equations llSj] are 



At = A 



'^Axx - if A, 
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(5) 

/t - fxx - \A\-j,. (6) 

An immediate prediction of (|4]) with the MC exponents, 
that the amplitude of solutions u{x,t) of ^ (computed 
as the time-averaged root-mean-square (rms) amplitude 
((u^)^/^), where (■) and ~ denote time and space aver- 
ages, respectively) should scale as ((u'^Y^^) = 0{e^/^) = 
0(r^/^), was verified in |ll,[li|; see [l^ for an analogous 
result in two space dimensions. 

We have performed the multiple-scale analysis of ^ 
to the next order in e, using the Ansatz (j4]) with the 
MC scaling a=3/2, /3=2, to obtain the asymptotically 
consistent 0{e) amplitude equations 
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(7) 
(8) 



(the linear 0{e) terms in ([7]) correspond to next-order 
corrections to the exact dispersion relation o'(fc)). Ob- 
serve that unlike the 0{1) equations (l5])-(|6]), these 0{e) 
equations also capture the roll solutions Uq {x) ^ of the 
GL regime, which have size \A\ = 0{e~^^'^) in this scal- 
ing. One would expect, though, that in the chaotic MC 
regime (in which u = 0{e^/^), \A\ = 0(1)), the 0{e) 
correction terms in (I7|)-(l8]) would lead to small 0{e) per- 
turbations to the solutions and their scaling. 



III. MODULATION EQUATIONS 

In the light of the decomposition ([3]) for solutions of 
(fT|), one may seek a modulation theory for the slow 
evolution of the envelopes, coupling the amplitude Ui 



IV. NUMERICAL COMPUTATION OF 
SCALING EXPONENTS 

Following previous studies [l^ [2l| , we set out to ver- 
ify that Ui and Uq indeed scale for small e as predicted 
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FIG. 4: (Color online) Matthews-Cox scaling of the pattern 
mode Ui, and anomalous scaling of the large-scale mode Uo'- 

(a) Scaling of averaged rms amplitudes of \Ui\ from ([1]) and 

|e^/2^| from (O-®; we find {{Julfy^^) ~ £^/2((|Ip)i/2) ~ 
for q' = 1.49±0.01. The (black) dashed line indicates the 
MC scaling prediction a = 3/2. 

(b) Scaling of averaged rms amplitudes of Uo from U]) and 

e\f from dZl)"®; we find {(cVy/^) - e''((72)i/2) ~ e"' for 
13' = 1.75 ± 0.03 (black dashed reference line). The (blue) 
dotted line indicates the MC prediction /3 = 2. 
Symbols for this and the next figures (hardly distinguishable 
here): (Magenta and red) Triangles: L — e£ = 2tt ■ 12.8; 
(blue) Circles and (black) Diamonds: L = — 271 ■ 25.6. >, 
o: Nikolaevskiy PDE £ = 2-n-m. ■ 12.8/e, m=l and 2; <, o: 
0(e) amplitude equations ([7l-([8]), L = 2nm ■ 12.8, m=l and 
2. 



by the MC exponents a=3/2, ^=2. With a view to iso- 
lating and thereby computing the individual scaling ex- 
ponents directly, Tanaka (2T| numerically integrated the 
PDE H]) on a domain of length ^=512 for g [0.001,0.4], 
and computed the rms scaling of the Fourier coefficients 
nearest in his system to k—0 and fc=l; these Fourier 
modes were used as proxies for Uq and Ui to argue that 
the results were consistent with MC scaling. However, 
the numerical evidence of [2lj, Fig. 5] does not appear 
sufficiently well-averaged to establish the asserted O(e^) 
large-scale scaling convincingly. 

Our implementation of Fourier filters to extract the 
full slowly-varying envelopes Uq and Ui from the com- 
puted solution u{x,t) of the Nikolaevskiy PDE ([T]) pro- 
vides a more powerful method to estimate the crit- 
ical exponents, via the e-dependence of the time- 
averaged rms magnitudes of Uq, Ui. In Fig. [3] we plot 

and ((J7o^)^^^) computed as functions of e 

^-5 



for e [10 ^,0.04]. As seen in Fig. Ufa), we obtain 

(([C/i 1^)1/2) = 0(e^) for a' w 1.5; that is, we find in 
agreement with |21| that the pattern mode Ui indeed 
satisfies the MC scaling a'=a=3/2 [T^. 

The computed scaling of the reconstructed large-scale 



FIG. 5: (Color online) Scaling of averaged rms amplitudes 
of £~^Uo from H]) and / from ©-((HI). This is the data of 
Fig. [IJb) rescaled vertically by e~'^ to emphasize the anoma- 
lous scaling: according to the prediction of the Matthews-Cox 
scaling Ansatz, which for the large-scale mode is / = 0(1), 
Uo = 0(2^), these values are expected to be e-independent. 
Symbols are as in Fig. 21 The solid black squares are obtained 
from additional computations of the Nikolaevskiy PDE Q for 
larger L ^ ^ 2tv ■ 51.2 (m=4). 



mode Uq is rather more surprising: as shown in Fig.mjb), 
over a wide range G [10^^,0.04] our data is consis- 
tent with {{Uo'^y^^) = Oie^^') with the anomalous expo- 
nent (3' w 1.75, contrary to the prediction f3'—f3—2 from 
multiple-scale analysis. In the alternative representation 
of our data in Fig. [5] the deviation from the MC predic- 
tion is more readily apparent: e~^((J7o^)^^^) is clearly not 
e-independent, as confirmed also by some larger-domain 
computations performed as an additional check. 

The leading-order amplitude equations ([S])-®, being 
e-independent, are manifestly unable to capture this un- 
expected scaling behavior. However, for comparison with 
the full Nikolaevskiy PDE, we have integrated the 0{e) 
amplitude equations for A and / as functions 

of the slow variables X G [0,L] and T, for the same 
e- and L- values as before. The functions u{x,t) recon- 
structed from A and / using ^ are qualitatively similar 
to computed solutions of ([1]); for e > trajectories ap- 
pear chaotic and are not observed to settle down. 

We extracted amplitude scaling behavior for (I7])-(|S]) by 
computing the time-averaged rms amplitudes of je'^^^^j 
and e^f; as shown in Fig. 21 the results agree remark- 
ably well with those obtained for \Ui\ and Uq from 
the full PDE (dl). Equivalently, the statistics of the 
pattern mode A are, unsurprisingly, asymptotically e- 
independent; but for the large-scale mode, we find that 
approximately {{PY^"^) = C'(e''/^), or, as shown in 
Fig. m {(f^y^^) = 0{e^^/'^): whereas one might have 
anticipated a small {0{e)) effect of the added 0{e) terms 
in the amplitude equations (I71)-(|S1), in fact they are able 
to capture the (asymptotically large) corrections to MC 
scaling. 

We observe in Fig. |5l a possible weak dependence of 
((72)1/2) (^j^^i ^-2(^^^)1/2')) ^jjg system size L = 
ei, which may merit further study; but that for each 
L we have the same scaling with e (with some possible 
finite-size effects appearing at the smaller domain size 
for e^ < 10~^). The deviation for e^ > 0.04, where 
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FIG. 6: (Color online) (a) Scaling of correlation times for the 
rms values of \Ui\ from ([l]), and \A\ from ©-((S]); we find 
n TA ~ 2.3, decreasing slightly for larger L. (b) Scaling 
of correlation times for the rms values of Uq from U]) and / 
from (O-®. Symbols are as in Fig. ID 



the separation of scales is weak and the "perturbation" 
parameter is relatively large, is unsurprising. 

The corrections to scaling appear not only in the rms 
amplitudes: While the multiple-scale Ansatz (U) intro- 
duces only one slow time scale T = e^i, a glance at time 
series of Fourier coefficients Uk (t) computed from (HI) for 
fcwO and fc«l, or at the evolution of Uq and \Ui\, suggests 
that for small e, at large scales the dynamics are much 
slower than for the pattern mode. 

To begin to quantify this observation of potentially dis- 
tinct time scales, for solutions of ([T]) we considered time 

series {Uo{-,Tfy/^ and {\Ui{-,T)fy/^ of the rms mode 
amplitudes depending on the slow time variable T, and 
computed (full-width at half-maximum) autocorrelation 
times tq, ti as functions of s. Again for comparison, we 
similarly computed correlation times tj, ta from time 
series of rms values of /, |^| computed via ([71)-([5]). Our 
results are shown in Fig. [6l again, the 0{e) amplitude 
equations closely match the behavior of the full PDE. 

The asymptotic validity of the leading-order modu- 
lation equations (l5])-(j6l) implies that these correlation 
times should be (asymptotically) e-independent. Indeed, 
Fig. inija) shows these predictions again satisfied for the 
pattern mode Ui: we find ti f» ~ 2.3, seemingly in- 
dependent of £ (though with, again, an apparent weak 



dependence on domain size L). However, as seen in 
Fig. MJ^) the long-wave mode Uq is substantially slower: 
while our statistics are not yet well- converged (particu- 
larly for small e), our data appear to indicate the pres- 
ence of another time scale diverging as e 0, and scaling 



possibly as tq w r/ 
may be near 4/3. 
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DISCUSSION 



The description that emerges for the Nikolaevskiy PDE 
([T|) with sufficiently small e is of spatiotemporally chaotic 
dynamics with strong scale separation and anomalous 
long-wave behavior. The solution space features two dis- 
tinct scaling regimes: The 0{e) GL regime, which con- 
tains the roll solutions Uq{x) together with their stable 
manifolds, is unstable to spatially varying perturbations, 
which lead to collapse to a smaller, apparently attract- 
ing, set in phase space, the 0{e^^^) MC regime of "Niko- 
laevskiy chaos" . 

Our observed anomalous scaling of the long- wave mode 
Uq implies that the leading-order amplitude equations 
do not fully describe the dynamics on the attrac- 
tor of Nikolaevskiy STC, at least for the parameter range 

G [10~^,0.04] we investigated. This conclusion is sup- 
ported by recent observations of non-extensive behavior 
in the Matthews-Cox equations (26| . 

The higher-order modulation equations, however, do 
appear to capture the corrections to amplitude and time 
scaling: the function f{X,T), evolving according to the 
coupled equations (IZl)-(|ll), scales in the same way as 
e^^Uo, where C/q is the large-scale mode extracted from 
solutions of (HI by Fourier filtering. Further investigation 
of (I7])-(|H]) is needed to clarify the mechanism by which 
the addition of 0{e) terms to the MC equations modifies 
the scaling. 
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